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ABSTRACT 



Tho Tcrzaghi Theory of Consolidation does not con- 
sidor tho tirao-dopondent loading conditions, oncountorod in 
ordinary construction practice, except in a graphical approx- 
imation, Tho development of a rigorous mathematical solution 
to this problem by Schiffman requires simplifying assumptions 
relative ro tho permeability and consolidation characteristics 
of the soil in order that the general differential equation for 
consolidation under such loadin: may be linearized, 

Tho investigation of tho validity of those assumpt- 
ions for the case of one -dimensional flow was conducted, 
utilizing both standard lover-arm and bellows type consolidation 
equipment. The assembly and use of a constant-head permoamotor 
to measure flov/ volume through a consolidating soil mass pro- 
vided, for tho first time, an opportunity to correlate instant- 
aneous and average values of tho coefficient of permeability, 
v/ith tho porosity oxlstant at tho time of reading. The wide 
variance in instantaneous values dictated tho use of finite 
incremental average values. By limiting tho time increment 
to 15 minutes, it was found that analysis of data was con- 
siderably improved. 

Pour consolidation tests v/oro conducted, using 
instantaneous standard incremental loading, timo-dopondont 
loading and small incremental loading, A proposed fifth test 

ix 



using small incremental loading r;as cancelled due to failure 
of a metal casting in the equipment under 16T3P loading. A 
pure Kaolin clay v/as tested in order to limit the number of 
variables entering into the problom. It v/as considered that 
the use of undisturbed natural soil samples would introduce 
intangible factors which would seriously interfere with the 
objectives of this pilot study. 

Analysis of tost data shows that, for construction- 
type (i.o. time -dependent loading), the assumption of linear 
variation of permeability with porosity, utilized by Schiffman 
in his solution of the general differential equation of con- 
solidation, is valid, but only for small load increments. The 
linearity of this relationship does not appear to hold under 
conditions of largo absolute loads applied over a short tliiio 
period. It is hypothisizod that the apparent exponential 
relationship is the result of the shearing of adsorbed water 
from the clay particle during an intense structural densi- 
flcation under such loading. Further study of this problom 
was strongly rocorainondod. 

The approximation technique utilized by Schiffman, 
in the solution of the variable loc rmoability case, was found 
to bo limited in accuracy only by the tine-increment used by 
the investigator, A wide variance in average values of the 
coefficient of permeability was found to exist under all 
conditions of loading, which factor requires almost constant 
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attendance if the average value is to approach the instantaneous 
value. 

No conclusion was drawn relative to the consolidat- 
ion-pormoability-timo relationship under time dependent load- 
ing, since fitting procedures for the tost curves wore not 
established during this investigation. 
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PART 1. 
li TRODUCTION 

A. Objective 

The furnishing of a reasonably accurate engineer- 
ing estimate of the settlement to be encountered in the 
founding of man's structures on soil strata which have been 
subjected to the vagaries and assaults of natural forces for 
raillenia of time, has been one of the most ir.iportant con- 
tributions of the field of soil mechanics to the Engineering 
Profession. The complexity of the solution to this problem 
cannot be overstated, since both interdependent and non- 
dependent, linear and non-linear, variables are encountered 
throughout any analysis of the settlement characteristics of 
Engineering Soils. 

It is the current practice of Soil liechanics and 
Foundation Engineers to compute total settlement and the 
time-rate of settlement of loaded compressible strata in 
accordance with the mathematical analysis developed by 
Dr. Karl Terzaghi in his theory of consolidation. Although 
rigidly correct for the assuiuptions used by Dr. Terzaghi, 
the use of the Terzaghi theory does not alv/ays render 
satisfactory results, since the simplifying assumptions of the 
theory cannot usually be found in the field. Among the 
assuraptions used by Dr. Terzaghi is the concept that loads 
are applied to the soil instantaneously v;lchout changing its 
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permeability characteristics. Although several others of 
the basic assumptions of the theory of consolidation are 
subject to critical revue, an extension to the Terzaghi 
theory as developed by Assistant Professor R,L. Schiffman 
of the Civil Engineering Department, Rensselaer Polytechnic 
Institute, proposes a solution to the settlement problem 
which considers the loading to be time-dependent and the 
permeability to be variable. Since this analysis closely 
approximates field conditions, it is considered that a more 
accurate appraisal of the action of soils \mder loading may 
be gained by its use, vice the original Terzaghi theory. 

With certain si.rplifying assumptions, the 
Schiffman concepts can be used to predict the total settle- 
ment and the time-rate of settlement of compressible strata 
under time-dependent loading, within the limits of 
mathematical accuracy and the validity of the siiiplifying 
assumptions. The validity of those assumptions is the 
standard against which the application of Professor 
Shiffman's solution must be measured, if the solution is to 
be considered a useful contribution to the practice of Civil 
Engineering, 

It is the specific objective of this thesis to 
conduct a series of tests designed to investigate the time 
effects of continuously varying loading on the consolidation 
of soils and the variation of soil permeability with pore 
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pressures induced by such loadinij. 

The general objective of theso investigations is 
to contribute to a more thorough and basic understanding of 
the action of compressible soils under loading, with 
commcnsuratcly improving accuracy in the prediction of 
settlements. 
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B. Historical Review 

The accurate prediction ox the settlement 
characteristics of compressible soil strata under loading has 
been a major topic for research and discussion for several 
decades. The statement of Collingwood (1) in 1891 that 
sound undisturbed earth "should be penetrated to a 
sufficient depth to insure that it is not underlain by semi- 
fluid or compressible material, which may in time yield and 
cause trouble and danger..." while qualitative in nature, 
nevertheless Indicates an awareness by practicing Civil 
Engineers of that ora, of the basic cause of the settlement 
phenomenon. 

Although the scientific approach to the under- 
standing of soil action had been in the process of develop- 
ment by such pioneers as Wollny 1879-98 (2), Schlictor 1397 
(3), King 1899 (4), et al it is to bo noted that their v/ork, 
v/hilo contributory to the broad base of fundamental soil 
knowledge, v/as primarily directed toward the agricultural 
uses of soils. Thus the use of engineering judgement and 
experience, and ornperically derived formulae remained the 
tools of the Foundation Engineer until well into the second 
decade of the present century. 

V/ith the publication in 1925 of his "Theory of 
Consolidation", Dr. Karl Torzaghi (5) provided a raajor 
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contribution to the scientific evaluation of one of the 
Foundation ilhginoors most perplexing problems - the predict- 
ion of settlements to bo expected after loading a compressible 
soil. Using fundamental physical lav/s and simplifying 
assuriiptions, Dr. Terzaghi developed a mathematical treatment 
of the coiiplcx functional inter-relationships of various soil 
properties which by thermodynamic analogue yielded a 
differential equation solvable by use of Fourier Series, 

His theory, based in part on instantaneously applied load- 
ing and constant permeability, for the first time provided a 
quantitative, albeit not completely accurate, measure of the 
consolidation of co' pressible soil strata under loading. By 
his solution of the Consolidation oroblom Dr. Terzaghi, in 
the words of Professor P.F . Tsehebotarioff, "became the found- 
er of the now science of soil mechanics" (6), In 1941 
Professor E,J. Kilkawley (7) brought together the fundamental 
concepts underlying Dr. Terzaghi 's work, and presented a 
solution of the differential equation under conditions of 
loading and drainage encountered in the field. 

It is obvious that the assumption of instant- 
antaneously anplied loading, vfnile achievable in the 
laboratory, cannot be duplicated in the erection of a 
structure. In an effort to account for this incongruity in 
computing time-rates and total settlement, Dr. Terzaghi 
advanced an approximate graphical method which after 
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extensions by Professor Gilbuy and Professor Taylor (8) has 
been accepted as the closest available approach to field 
loading conditions. The recent work of Professor Schiffman 
(9) in developing a rigorous raathoraatical analysis of the 
consolidation process, as it is affected by the time-rate of 
loading, permits an analytical evaluation of the settlement 
to bo encountered without recourse to graphical approximations 
It is the purpose of this thesis to investigate the validity 
of Professor Schiff ..lan' s assumptions in extending Dr, Torzashi 
v/ork. 
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PART 11. 

TH30RY 

A. The Consolidation Process 

The compression of soils under externally applied 
loads, either natural or nian-inade, is the basic cause of 
sottloncnt of structures. This compression is intimately 
associated with the pore spaces, (ie. intergranular voids) 
of the soil structure, and the escape of ground water, there- 
from, Since both water and the soil grains are considered in- 
compressible, the only way that the height of a soil mass can 
be reduced, in the vortical plane, is by the escape of the 
pore v/ator from the soil structure v/ith a concurrent, 
donsifying, structural rearrangement of the soil grains. 

Thus it is scon that the rate of compression of a soil is a 
function of tho rate of cscanc of pore water. 

With the application of load to a saturated 
perious soil mass such as clean sand, the escape of pore water 
is almost instantaneous, since the pcrvlousnose of tho soil 
places no obstacle to free passage of water. On the other hand, 
if tho loaded soil is a saturated clay, the escape of tho pore 
water requires considerable time for completion. At the 
instant after load application, before any pore water has es- 
caped, the soil structuro cannot support any of the applied 
load, since it is in equilibrium, and any change in this con- 
dition requires rearrangement (io. compression) . It follov/s 



8 



that the load is supported, toinporarlly by the pore \vator, 
which must incur an incrcaso in prossuro to perform this 
task. This increase in pore water pressure is termed 
"hydrostatic excess pressure". Under conditions of boundary 
draina£;c, the hydrostatic excess pressure is subsequently 
relieved by the escape of pore water. This pressure relief 
by boundary drainage, requires transfer of stress from the 
escaping pore-water to the grains of the soil structure, which 
reacts to the induced forces of loading by a reduction in 
volume oquivllcnt to the volume of escaping pore water, thus 
permitting additional intergranular contact to absorb the 
transferred stresses. In this manner, the stresses duo to 
loading slov/ly pass from hydrostatic excess pressures to 
intergranular or effective pressures, v/ith full stress assujiipt- 
ion by the soil structure, upon complete dissipation of the 
hydrostatic excess. This adjustment process whether it be 
relatively fast, as in sands, or very slov/, as in impervious 
clays, is termed "consolidation". In its simplest form tho 
consolidation process is reducible at any instant of time to 
tho equation; 

p -- u t- 0 " 

where p is the comprossivo stress duo to the applied load 
u is tho hydrostatic excess pressure 
Z~ is tho offoctivo pressure 

For tho limiting conditions of (1) no dissipation of hydro- 
static excess and (2) complete dissipation of hydrostatic 
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equation reduces to 



(1) 


P.- ^ 


(2) 


P-=.'7, 
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3. Tho Torzap;hl Theory of Consolidation 

The analysis of the consolidation process develop- 
ed by Terzaghi ia based on the iollo\7ing assur.iptions ; 

1, Homogonoous soil mass 

2, Complete saturation 

3, Negligible compression of soil grains and v/ater 

4, Action of inf initossimal masses no different 
from that of larger, representative iaasses 

5, Onc-dimcnsional compression 

6, One -dimensional flow 

7, Darcy's Lav/ is strictly valid 

8, Constant values for certain soil properties 
which actually vary with pressure 

9, Void ratio varies linearly \/ith apnlicd 
pressures 

»/ith the exception of the ninch listed assumption only minor 
inaccuracies are considered bo stoi.i from these assumptions. 

The linear relationship between void-ratio and pressure in- 
troduces serious discropencics 'diich limit tho valadity of the 
solution. However, the solution would bo \.'Cll-nigh impossibly 
complex under any other type of fu*ictional variance. 

Figure 1 illustrates the consolidation process under 
consideration, A clay layer of thickness 2H, between two 
pervious sand layers is stressed by an applied unie load p. 
Boundary Drainage permits dissipation of tho hydrostatic 
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excess pressures induced in the pore ’./ator. As dissipation 
progresses the effective pressure & increases from 0 at tiino 
to: p at tinio when the process is completed. At inter- 
mediate times -etc . the consolidation process is shov/n to 

be more advanced tov/ard the surface and loss advanced at the 
center of the stratum. Consider, now, a differential prism 
of soil from the upner half of the clay layer, with cross 
sectional area of 1 and height dz as shown in Figure 2. The 
piezometer tubes are considered to accurately measure pore- 
water pressure at top and bottom boundaries. Since pore 
water is escaping through the upper boundary, a drop in head 
is occuring in the direction of flow. This drop in head, dh 
is related to the dissipation of hydroscatic excess du by 

dh -r ^ vh :runit wt of water 

Ow “ 

the hydraulic gradient, i, is defined as the drop of head over 



a given distance. 



dh 



By substitution 



du 



Darcy’s Law for flow of water through soil masses states that 

the rate of flow is proportional to the hydraulic gradient 

V - ki k r coefficient of 

permeability 



V - flow velocity 
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Again by substitution 



V _ k du 

By differentiation, the change in flow velocity over distance 
dZ during a given time interval dt equals 

^ _ U^yL^. 

Since Darcy’s Law can also be stated in the form: 

V - Q where Q - discharge 

At 

A = Cross Sectional Area 
t = time 

it follows that, in time interval dt v/ith A- 1, V-Q (i.c.: that 
the flow velocity represents the amount of water entering the 
base of the prism sketched in figure 2) then, the increment 
of velocity gained over the height of the prism, dz, must equal 
the increase in discharge through the top of the prism, ic. 



dv = dQ 



Any increase in discharge from a fully saturated 
soil must result in a decrease in pore volume. This Porsity 
change .6n can be expressed as; 



An ^ -A^-. 

/ *e 

and since the coefficient 



whore e - void ratio 



of compressibility a^ 



AP 



/An = 

e 



t 
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\i'hcn the change in pore volui'.io A.'i is completed ^ the pressure, 
p, equals the effective pressure carried by the soil grains, ''T, 
and the preceding equation can be expressed as 

_ d o' 

dt di 

Any increase in Q under constant unit load P during time 
interval dt must equal the decrease in Hydrostatic itccoss 



Pressure, U, during the same 


time period* 




^ u, da 




3n _ 


dm dv 


ir ■ Ji' 




di 


dt !+e 


now since 








A-- - - - -Ijt- 






- - ^ 






hi 


di 


and 

an 


d\X 






ht 


U 


!*e 




it follows that; 
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"at" ' 


The fundamontial differential 
equation of consolidation 


where Cv . k(i + c ) - 


= Coefficient 


of Consolidation 











0 " 

The solution of the funda;:icntal differential equation by 
Fourier Series, for the boundary conditions as shown in 
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Figure 1, is as follows 






iLll''* 

A \ 



n TO 



^ tTo 

Whero a dimensionless Time factor 



5 iK 



T 



lZn*i)7rz 

Z H 

- t 
H 



(2aM)^77^T 
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With Cv - Coefficient of Consolidation 

t ^ time for Hydrosto.tic Excess to decrease to U 
H •= Longest Drainage rath for v;ator flow 



By plotting values of U for various values of and T ^ 
the family of curves shovm in Figure 3 eliminates the tedious 
solution of the consolidation equation and permits the use of 
Figure 4 as a practical method of arriving at a prediction of 
settlement under loading. Since the ultimate settlement 

- iU Ap 

I t-e 

then the settlement at any time interval equals 






V/hore Uavg =■ Average Hydrostatic Excess 
Ui - Initial Hydrostatic Excess 



Hov/ever, among the factors ignored by the Terzaghi 
solution is the gradual > rather than instantaneously > applied 
loading encountered in the construction process. A graphical 
approximation, based on the premise that, at the end of con- 
struction, the settlement is the same as that which would have 
resulted in half the time, if the entire load had boon 
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applied throughout, v/as first proposed by Torzaghi and Gilboy. 
Taylor extended it to provide for predicting sottloinonts dur- 
ing the loading period, by assuraing that at any specified 
percentage of the loading period, the load acting equals this 
porcontago of the total load; and at this time, the sottlomont 
equals the sottlomont at ono-half tho time in question from 
the curve of time vs, sottlomont for instantaneous loading. 
Figure 5 illustrates this correction. 

The solution of tho fundamental differential equat- 
ion wfhich by tabulation gives tho curves of Figures 3 and 4, 
was based upon honogonoous boundary conditions, i,o,, 
prossuros at the top and bottom of tho sample are equal. Tho 
use of a constant head perracametor imposos different boundary 
conditions than these, in that tho applied load is rosistod by 
the constant pressure-head of tho permoamoter, Tho following 
solution to this problem has boon developed by Professor 
Schiffman and graphical values for tho loads applied are pre- 
sented as Figures 6, 7a and 7b: 



v/hcro 
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u - 
u. 
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1 - Consolidation 

Pressure duo to Pormoaraoter Head 
pressure duo to applied load increment 
Average Hydrostatic Escess Pressure 

^ Consolidation with Homogonoous Boundary Pressures 
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As may bo noted from Figures 6,7a and 7b the solution of the 
equation gives the same valuo for the Tima Factor as that 
tabulated by Taylor (8). The basic assumptions of linear 
variation of hydrostatic excess pressure with depth thus 
reduces a seciaingly involved problem to the use of tabulated 
values, instead of tedious computations for each variation 



in Permoameter head. 
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C . Schlffuan Sxtonsion to the Theory of Consolidation 



Professor Schiffnian* s development of a three- 



dimensional anisotropic theory of consolidation ( ) provides 

a working basis for the study herein conducted. His assxunpt- 
ions in developing this theory are; 

(1) Complete saturation of Soil Mass 

(2) Incompressible fluid 

(3) Incompressible soil solids of small particle size 

(4) Darcy’s Law valid at any instant of time 

(5) Change in volume is linear with imposed pressure 

(6) Change in voliime is small compared to original volume 



S, , In accordance with the ]a w of conservation of mass, the 
fluid flow into the r,iass through the surface must equal the 
total volume change of the mass. In terras of vector functions 
the total instantaneous mass flow is equal to 



Consider a prism of soil of volume, V, and surface, 




where 



j V’l'ldS “ flow into mass through surface S 



and 



/ 




increase in volume due te 



internal flow generation at 
rate, Q, due to head 
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the voliuno change duo to flov/ is equal to: 



()t 



\ 



0 V./ ^ \/^ 

■ 5 T 



dv = o 



Equating, v/o arrive at: 

j V - ( hx- h) > Q ^ 

SincG the voluiTD is arbitrarily cstablishod vi?c can v/rito 

V'(kvh)tQ -- 0 

Since only tliosc strains duo exclusively to volume change are 
herein considered (assumptions 5 and 6), only normal stresses 
need be treated. Tho i[odulus of Volume Change, is defined 
by: 



hrv = 



de 



V/hcro 



e = (r. . (T, . or. 



Since the total pressure or the Soil Hass is composed of the 
neutral pressure, u, and tho effective moan stress 6 we can 
wr i t G : 

6 ' 0 ^ UL 

/here u is the pore pressure defined by: 

VO - 

Since che consolidation process involves the 
transfer of initially incurred stress from the water to the 
soil grains, we can v/rito: 



6u. 
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and that: 



()vy - * tn ^ u. 



By substitution in the equation of mass conservation, 
v/e then arrive at the differential equation of consolidation 
for a variable permeability and time-dependent loadirg : 

V • 7U.) 1 - q yio = m 

the treatment of the two dependent variables, k and u, is based 
upon the work of Schmid { ), v/ho proposes the relationship: 

k /S (n-n^) 

whore /3 s Constant 

n = Porosity 

•= Ineffective Porosity 



Using this relationship, wo can extend it to say that a linear 
relationship also exists between k and u: 



v/horo 



k 



- c>c ec t 





z Coefficient of Consolidation at the end 
of consolidation 



and 



(X r Modulus of Permeability Varation 
Oe 5 V O'! 



It is now feasible to establish a differential 
equation covering the general theory of consolidation under 
conditions of varying permeability and time dependent loading; 




I 
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(a) In terms of excess pore pressure u 




(b) In terms of the permeability k 




( c ) and 



-<X.(V<A«- 7t(.) + CiiCVu>)*- .(h:^-L(XL(.o')VV 



The present study being limited to the considerat- 



ion of one-dimensional consolidation, the general equation 
for such a case can be written; 



R -r Rato of change of imposed excess pore pressure 
C ■= Coefficient of Consolidation at start of con- 
solidation. 



This equation can be written under the following conditions 



(1) Double Drainage 

(2) Infinite extent of soil mass horizontally 
(5) Finite thickness of soil mass 

(4) Uniform initial Coefficient of Permeability 

(5) Uniform initially imposed pore pressure 





whore 
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Professor Schiffman has prosonted several individual solutions 
involving the following problems; 

(1) Constant permeability and general time dependent loading 

(2) Constant permeability and linear loading 

(3) Constant Permeability and construction loading 

(4) Constant Permeability and harmonic loading 

(5) Variable permeability 



His solutions to these problems are as follows: 

(1 ) Constant permeability - time dependent loadin g 

u(Ht) = -i- mz) njii A)££ 

^ ^ J.H “I 

* IT ^ \i r (^ Vdv 

ni L'oU:, ZH J 



( 2 ) Constant Pornoabilty » lire ar loadirg 

_i_. imonrr/. ^-sdiLlr) 

ToTr3 L ni Tu ^ f 



ri iJiT- 

Since the average pore pressure can be expressed as; 

■/o 

Then by substitution 

U(T): Mis C -J_ ^ ) 

Whore T is a dimensionless Time Factor 

T =-- t 

Figure 8a shows the curve developed from a computation 
based on this equation, from which all other computations for 



the linear case can be made 



za 




TIME FACTOR (T) 
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riguru 8 (b) proscntij the computation involving the 
linear loading period ending at time t and the amount of 
excess pore pressure dissipated at that time. The oercontage 
of settlement completed is oquivilont to the amount of oxc>^ss 
pore pressure dissipated and thus the settlement at the end of 
a linear type construction loading period can be computed. 

As an example of the usefulness of Figure 8b, let 
us consider a structure founded on a clay structum VTith Cy'l-5^ / 
If, then, the length of construction is related to the end 
time factor To; 

by entering Figure 8b, a table of percentage cf settlement 
dissipated can be constructed, in terns of the length of 
construction. Table 1 shows this construction as developed 
by Professor Schiffman, 

By use of such tabl.s, the engineer can now ra ke a 
rational decision as to the desired construction time period 
which will restrict post-construction settlements to satisfact- 
ory values, or as to solv^ction of alternate solutions, 

( 3 ) Constant Permeability-Construction Loading 

With the results of die previous sections available, 
the pore pressure can nov; bo analyzed at any instant and for 
any construction period. Figure Be presencs a typical load- 
tlme diagram for construction loading, wherein a load I, is 
imposed at time t, , the end of a linear construction period. 
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Treatment of this problem is as doscribod in (2) above. 

Figure 8d and 8c present the solution of the 
determination of average pore pressure during and after con- 
struction, These curves can be used to determine the 
theoretical consolidation curve in the following manner: 

(1) Estimate the construction tine 

( 2 ) Determine the Coefficient of Consolidation from 
Laboratory tests 

(3) Dctcrmi:co To from the Coefficient of Consolidation 

(4) Enter Figures Od :aid 8e and select proper Con- 
solidation Curve 



( 4 ) Variable Formoabil 1 b:/ 

The basic equation governing thiS case is; 



Ic \ SLi^] 4 ^ - 

na 



JSS.. 



U. 









m 









However, the solution of this equation is hariiporod 
by its non-linearity, for '-rhich general techniques of solution 
arc unavailable. Solution bj approximation techniques is 
thcroforc utilized, in lieu cf numerical computation by 
analogue or digital computer methods. 

The incremental time approximation is first con- 
sidered, It is assumed the per'icability remains constant 
over a finite time increment, but that the permeability varies 
from one incro'.ront to the next. For a finite number of time 
increments we have: 

ua.r'i- ^ Y_ 

1^* - 



1 



3V 



1 . 0 , 

0.9 

0.8 

0.7 

06 

05 

0.4 

0.3 

0.2 

0.1 

0 




TIME FACTOR (T) 






■t 





3S 




TIME FACTOR (T) 



U.( T') 







36 



r‘ - 



/ 






By fitting incroi'cntal laboratory tost curves to the 
theoretical curve, on the basis of theoretical values of u 
and t', successive values of the bracketed term can be 
detornined. The degree of accuracy desired can be controlled 
by the number of incroi.ionts utilized in the fitting procedure, 
A second approximation considers that an ex- 
ponential permeability pore-pressure rolationship will 
approximate the linear condition originally assiU'.aed, over the 
entire range of the consolidation process. The governing 
differential equation is: 

6 






7ith boundary conditions: 

(-) 

(g) ktx.o'k 






- auo + 



The solution is: ^ g ^ 



which when converted back to excess pore pressure is: 

mz.t -L s,„ JXB bt 

I.3i5’ - 

where B is the Coefficient of Consolidation Pernoability and 



is c qual t e : 
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where 

h ' U ^ 

The solution to the general case for ono-dimonsional 
consolidation can bo utilized for varying permeability, if a 
now Time Factor V is substituted for T whore: 



V ~ -fi- J 



It is noted that as approaches the value of B 

approaches Cc 

It is essential, in this solution, to satisfy the 
condition thab the soil have an initial and final coefficient 
of permeability, v/hich is constant throughout the mass and at 
terminal points in any segment# 
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D. Poriiioabllity 

Tho property of soil, concornud with tho facility 
of travel of water through a soil nass, is tormed Pornicability , 
Since the physical properties of soils as well as tho state of 
stress of tho pore water have important effects on the 
mechanical behavior of a soil mass, it can bo seen that this 
property is of major interest in tho study of consolidation 
characteristics of soils. 

The extremely small pore sizes encountered in most 
soils cause any flow of water through the soil mass to be 
laminar in nature, (i.o. the a.iount of head lost in friction 
is directly proportional to tho velocity of flow). Such 
laminar flow takes place in accordance v/ith dynamical 
eqatiens of motion which by analysis of various simple 
channel shapes have been reduced to usable formulas. In 
engineering problems the flow through individual flow channels 
is not required, but it is rather tho average flow through 
tho soil mass which is important, 

H, Darcy in 1856 demonstrated experimentally that 
the velocity of poro-v/atcr flow through a soil mass is directly 
proportional to the hydraulic gradient thus formulating tho 
basic law of flov7 

~ *• whore v - Velocity of flov/ 

k vr Coefficient of 
Permeability 

i Hydraulic Gradient 



I 
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Tho flow velocity, thus defined, is the superficial velocity 
through the soil mass, not tho specific velocity through each 
flov/ channel. 

To evaluate the coefficient of Permeability 
directly, the use of a constant head permcameter provides a 
quick, siriiplo solution. Consider a soil sample with cross- 
sectional area. A, and height, L, Let there be a supply of 
v/ator at constant elevation, II, above tho sample, and let there 
bo some moans of measuring flow volume through tho sample. 

With these values knovm. Darcy ' s Law can be re-written: 

V - q = ki whence ^ ^ 

A lA AH ■ AtH 

where Q = Total Volume of Plow and t = time, 

Tho use of a constant-huad peri-camoter has boon 
restricted to soil typos of relative coarse grading, since 
flov/ volume through fine grained soil is small. However by 
continuous moasuromunt of flow volume during a consolidation 
test with che apparatus shown in Figure 9, flow volume 
measurement is a simple matter. In addition to the obvious 
advantages of consolidation- time values for flow volume, this 
apparatus eliminates the problems of air-pressure application 
encountered by Peterson (12) in using a variable head permea- 
moter to measure flow through an Illito sar.mlo. In addition 
it makes possible, an almost Instantaneous evaluation of 
permeability at any small finite time liiterval. 
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Tho evaluation of the coefficient of permeability 
obtained by use of the apparatus described above is a function 
of (1) tho porosity of tho soil mass (2) tho shape and size of 
tho voids and (3) the density and viscosity of the fluid. 

Since we aro dealing with a pure Kaolin Clay, laminar flov/ may 
be assumed. In dealing with tho porosity of fine grained soils, 
it must bo noted that tho flov/ channel cross-sectional area is 
reduced by tho amount of ionically '‘bound" water v/hich is 
adsorbed on tho surface of tho soil particles and, being im- 
movable, blocks part of tho flov/ path, 

Tho work of Schmid (13) indicates that the effect 
of this "bound" water can be reflected in an extension to tho 
well-known Hagin-Poiscullo equation for laminar flow through a 
tube of constant cross-section (8), Schmid's equation; 

Vi/hcro y^runit v/t, of water 

Dj. = Effective Diameter of Soil Particles 
Vl r Porosity of Soil 

Y\„ ST ineffective Porosity i,o, channels blocked 
by "bound" water 

^ coefficient of viscosity of water 
indicates that tho coefficient of permeability is directly 
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proportional to the ''effective'’ porosity, (n - no). This 

concept would appear at variance v/ith the experimental data 

compiled by numerous investigators, which shows that plot of 

k vs log o as a straight line, i.e. an exponential 

relationship, Hov/ever as Schmid points out, the basic 

relationship, n ^_e_ is actually only an approximation of 

HC 

the series expansion of the function log 1 or log oh- 1 , which 

n e 

ignores the higher order terms. In actuality the relationship 
between n and e is logarithmic and thus Schmid's work can 
be reconciled with other investigators’ work if this factor 
is recognized. The error involved in ignoring the higher 
order terms is less than lO/o and decreases rapidly with in- 
creasing void-ratio, 

Tho straight-lino relationship between porosity 
and the coefficienb of permeability does not hold in the 
region ii'sjno and Sclmnid indicates that such deviation can 
be expected, since in deriving his equation he made two 
approximating assujnptionsj 

(1) That tho effective diameter, Dg , is constant. 

(2) That the ineffective porosity, n^ , is constant, 
which arc never achieved in actual testing. It is further 
suggested that the consolidation process, by virtue of its 
reduction in rate of porosity change at largo loading 
increments, v/ould indicate that such deviation should be 
expected in such regions. Since these regions are tho same 
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as where n -- no, The 3trai;;_;ht-line relationship would thus 
bo invalid for larger load incror.aonts wh^ro porosity changes 
are small relative to the applied load, 

The accurate evaluation of the coefficient of 
permeability, during the consolidation process, is of oxtromc 
importance, since the simplifying assumptions of both 
Terzaghi and Schiffnan, include k as a constant for the one- 
dimensional drainage conditions considered herein. The 
determination of the nanrior of variance of k v/ith decreasing 
porosity under consolidation loading would be of assistance 
in the solution of the moru complicated timo-and- space 
dependent consideration of the consolidation process as pro- 
posed by Schiffnian. Toward this end the evaluation of k 
during the consolidation process will bo attempted as part 
of this work. 
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PA^T 111 

MATEiHALS AND APPARATUS 

The soil samples used in the conduct of this 
investigation v/ere prepared from oven-dried pure Kaolin clay 
secured from vVARD'o NATURAL SCIENCE CORP., Rochester, N.Y. 
under the sample designation '‘Kaolinite-Dry Branch Ga. 

Dana >1^92 ", This relatively pure clay mineral exhibited a 
very sensitive reaction to moisture, as noted during the 
tests for determination of the liquid limit, where the addition 
of a few drops of distilled v/ater changed the blow count by 10, 

For the conduct of the control test, Number 1, a 
standard Fixed Ring Consolidometer as shown in Figure 10 was 
utilized. 

The application of time -dependent loading to the 
samples during tests Number 2, 3, 4 and was accomplished 
by means of the apparatus shown in Figure 11. A Conbel Model 
No, 350 Consolidometer was adapted for use v/ith a King 
Manufacturing Company visual bleed. The visual bleed, re- 
ferred to throughout as a time-dependent loading device, was 
filled to within approximately one and one-half (l-o-) inches 
of the air outlet v/ith a high-quality brake fluid. The 
pressure accumulator, and bellows on the Conbel equipment v/ere 
also filled with brake fluid in accordance with the manufactur- 



ers instructions. 
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For the direct measurement of permeability during 
all tests conducted, a constant head permearaeter as shown in 
Figure 9 was constructed and utilized. 
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PART IV 

Testing; Procedure 

A. Preparation of Sample 

All samples were prepared at the Liquid Limit de- 
termined in accordance with AoTil "Standard Method of Test for 
Liquid Limit of Soils", ASTM Designation D-423-39, Plow Curves 
used in determining the Liquid Limit are shown in Figure 12, 
Reproducibility of the Liquid Limit in mixing the sample prior 
to testing was considered to be satisfactorily attainable by 
use of the Standardized Spatula ilethod, as may be evidenced 
by the moisture content deviation being less than for all 
samples. The electric mixing of samples as proposed by 
Edmonds and li/arren (10) was attempted but not utilized since 
the entrainment of air v/ithin the sample appeared excessive. 

The sample v/as placed as follows: a 200 gram, 
dry weight, sample v/as prepared at 3;5 above the Liquid Limit 
and transferred from the evaporating dish to an open plastic 
cylinder 4 inches in height and 2,500 inches in diameter. The 
cylinder was sealed by a metal plate at one end and beaten 
strongly against a solid surface, concurrently rotating and 
tapping the sides of the cylinder with a metal rod. Entrapped 
air is visibly removed during this process. The saiTple is then 
extruded by means of a porous stone piston on to a glass plate 
covered with a filter paper. The planar surface, produced by 
extrusion into the glass plate, is then forced into the fixed 
ring and the sample accurately trimmed to the ring size. 
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B. Standard Consolidation Test 

Test No, 1 v/as conducted on a standard Fixed Ring 
consolidometer as described in Part 111, As suggested by 
Professor Burrnister (11) this tost was used as a pilot study 
to establish techniques of sample preparation and testing. 

Since no apparent swelling resulted from Immersion of the 
sample, and the strain did not exceed ,02 in, per in, it was 
decided that the testing cycle shown in the rav; data sheets of 
Appendix 1 would be satisfactory. It is to be noted that there 
is a dearth of literature on the preparation of disturbed 
samples for consolidation testxng, which factor makes a pilot 
test of this nature an absolute necessity for any testing of 
such materials. Appendix i also includes raw data from a 
concurrant pilot test conducted by V. I.IcGuffey on the same 
material with the same equipment, McGuffoy's results show 
that reproducible results can be attained if technique and 
procedure are carefully analyzed and correlated. 
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C . Tine Dependent Loading; Tests 

Tests No, 2 and 3 v/cre conducted usinj the apparatus 
described in Part 111, The Conbel equipment was calibrated 
by the use of dead weights vs. air pressure and it was determin- 
ed that the calibration charts furnished by the manufacturer 
v/ere accurate to within one pound of load applied to the 
sample in the 1-2 ton range. It is suggested that future 
Calibrations of this equipment be accomplished by a proving 
ring v;ith a planar upper surface to insure axial loading. 
Attempts to calibrate the equipment with standard laboratory 
proving rings were iinsuccessful, due to the tilting of the 
bellows at pressures above 4 psi due to lack of axial re- 
sistance to load by the proving ring. 

Figure 13 shows the calibration procedure followed 
for the time-dependent loading device. Tables 2 and 3 in- 
dicates calibration pressures and data in tabular form. 

Figure 14 indicates the operational procedure 
followed in applying time dependent loading to the sample 
by means of the tij.ie-dependent loading device. 

Raw data showing pressure increments and loading 
cycle times are contained in .Appendix 2 and 3, 
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D. Pseudo Time Dependent Loading Tests 

Tests 4 a ' Rd -5- v/ere conducted in an effort to 
approximate time -dependent application of load by use of 
small incremental loadings over small finite time intervals* 
Load v/as applied instantaneously by by-passing the time 
dependent loading device as shovm in Figure 11, Raw data 
showing pressure increments and loading cycle tiities are 
contained in Appendix 4 and "6 . 
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E, Permeability Measurements 

Direct measurement of Permeability was achieved 
by means of a constant head permeater as shown in Figure 9. 

A capillary tube was calibrated to ,1 c.c. and voluiae decrease 
noted at each reading of the consolidometer . Thus a direct 
relationship between permeability and void ratio was 
established for all points of the consolidation process. 

Considerable difficulty was experienced in 
establishing valid permeability data under small loading 
conditions (i.e. t.s.f,). The per.ieameter head is sufficient 

to cause sample uplift prior to application of load and as is 
shown in Appendix 4, affects the consolidation process to a 
considerable extent. It was therefore necessary to restrict 
permeability measurements to loadings of greater than -§■ t.s.f. 

It must also be noted that the direct measurement 
of k for such a material as Kaolinite involves extremely 
small measurements of flow volume. Such measurement is 
feasible under time -dependent loading conditions since the 
porosity is decreasing at a rate which permits accurate gage 
readings for flow volumes of sufficient magnitude as will 
minimize interpolation errors. However in instantaneously 
loaded samples, the load application: 

(1) Creates an initial surge of pressure which 
causes a negative flow into pormoaraoter for 
a period of approximately 3 minutes. 
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(2) 


Seriously restricts the evaluation of 
permeability data after the initial 3 
minute surge period, because of the lack 
of kno' ledge as to the time or pressure 
limits of its influence after the initial 
3 minutes. 


The problem 


involved in (2) above is reflected in widely 



scattered permeability values for all load increments of 
Tost IJo. 1, 
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PART V, 

RESULTS and DISCUSSION 

A. EQUIPMENT 

The test equipment, as designed and assembled, 
provided data of sufficient accuracy for a pilot study of 
this nature. For future development of tost data, requiring 
reproducible results for purposes of definitive analysis 
and for publication, the equipment should be modified to 
provide the following: 

(1) Precise control of load 
application during the Time 
Dependent Loading cycles, 

(2) Smaller bore, oquivilent volume 
permeameter to provide more 
precise calibration of flow 
volume . 

DISCUSSION: 

(1) The application of Load at a fixed linear 

rate is of extreme importance in the conduct of 
tests designed to investigate the validity of 
Schiffman's extension of the theory for the one 
dimensional flow condition, It is noted that the 
rate of change of imposed excess pore pressure, R, 
found in the governing differential equation, is in 
fact the rate of application of load. Therefore for 
any of the cases considered in his treatment of the 
one-dimensional flow condition, the loading device 
must be capable of precise control, if correlative 
data is to be developed. 
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(2) The use of the Constant-Head Peri'iieamotor is 

v/ell adapted for precise measurewont of permeability 
during the consolidation process. However, the small 
flow Volume, encountered with the use of clay soils, 
dictates a finer calibration of the horizontal tube, 
in order to minimize approximating errors. A smaller 
bore tube would provide the precision required but 
has the disadvantage of reduced volume availability. 
Such reduction in volume can only be compensated for 
in two ways: (1) Extreme length of tube or (2) Con- 
stant attendance to preclude evacuation of the tube. 

A reasonable balancing of these two factors would 
appear to be feasible. The increase in accuracy, 
resultant, cannot be estimated, but v\rould be of ex- 
treme value in analyzing permeability variation. 
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B, TEST NO.l - STANDARD CQNSOLIDATIOM TEST 

Test Data in graphical form is presented in 

Figure 16, 

Discussion; 

The data developed is typical of the "type'' curves 
for the Kaolin clay usedi The work of previous investigators 
(10) closely parallels the results achievedi The variance in 
values is reasonable and the shape of curves is similar. 
The problem of sample uplift, due to the head of 
the permeameter, precludes permeability measurements prior 
to the application of -g- ton per square foot loading. This, 

however, is not a serious difficulty since it should be noted 

) 

that the loading range involved is usually on the precompress- 
ion portion of the e-log P curve* Since usable data can be 
developed only on the straight-line portion of the curve, 
where the uplift problem is no longer cxistant, the problem 
does not actually affect test results. 

The determination of the permeability of the test 
samples, under instantaneous load application, is a problem 
which requires a great deal of interpretative analysis during 
the primary consolidation period. It is noted that an initial 
negative flow into the permeameter persists for a period of 
2-3 minutes, immediately follovs/ing the application of load. 
This flow can be interpreted in three ways: 
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(1) That the flow is due to the squeezing out of Pore 
water by the application of load, which process is 
essentially complete during the primary consolidat- 
ion period. V/ith the dissipation of such pore 
water, the head of the permoamotor then induces a 
positive flow. The validity of the initial positive 
flow measurements is subject to serious question 
since there is still existant a varying amount of 
hydrostatic excess, 

(2) That the flovif is due to a hydraulic balancing of the 
load-in^osed pore pressure and the permeameter head* 
Such balancing is subject to a time lag since the 
transmissibility of pore pressure differs at the 
top and bottom boundaries due to the counterhoad 

of the permoamotor, 

(3) A combination of (1) and (2) above. 
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C , TESTS MO. 2 and 5 - TII.C DEPEIID^NT LOADIIlG 

Test Data in graphical form is prosonted in 

Figure 17, 

DISCUSSION: 

The rate of application of load in tho majority 
of incremonts was such that tho consolidation process 
vms essentially concurrent v/ith load application. The 
curves developed show the same shape and range of 
Variation for each test. 

The increase in the permeability of tho sample, 
in the later stages of tho consolidation process, is a 
phenomenon which appears to offer a subject for future 
study. This increase is not predictable^ nor is it 
always present. It would appear that, at some stage of 
the plastic range, there is an internal adjustment of 
pore-channel '‘effective" areas, which causes scattered 
increases and decreases in the permeability. The 
effect is such that it has been termed the "yo-yo 
effect". The only reasonable hypotheses proposed at 
the present time are: 

(1) That dissolved gases in tho permeant 
(ie, v/ater in these tests) are condensing and provid- 
ing an increase in flow volmio, which is measurable 
on tho pormoamotor. 



(2) That tho ionically bound water "hull" which 
restricts tho flow of tho permeant during tho primary 
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consolidation poriod, is, by pressure adjustment, being 
sheared from the clay particles, thus increasing flow 
by the volume of the sheared hull and the increased 
volume of flow thru the enlarged pore-channel. De- 
crease in flov/ could be the result of structural re- 
arrangement, due to particles shifting, in response 
to the pressure void left by the shearing of the com- 
pressed water hull. The increased stress on the 
particles, v/hich would cause such shifting, is avail- 
able from the dissipation of the structural pscudo- 
hydrostatic-cxcoss pressure of the water hulli 

(3) That errors in guage readings have occured. 

Due to the cuinulativo type readings takon> and the 
repetition of this effect, it is considered that errors 
duo to (3) above are minimal. 
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D, TESTS NO. 4 - PSIHUDO TILG-PaPEND.^NT LOADING 

Test Data in graphical form is presented in Pigviro 18. 
Although it was hoped that another test of this nature could 
bo conducted, failure of a casting in the pressure lino from 
the bellows to the guago system precluded conducting it. 

DISCUSSION: 

The influence of the magnitude of the loading 
incromont is quite apparent in data developed by 
this test. The deviation from standard load in- 
cromonts does not unduly affect the c-log P curve 
but it is in the k-n relationship that the effect 
is most apparent and requires detailed analysis. 

As may bo soon from a study of Figure 17, the linear 
relationship botvi;oon the permeability and the porosity 
is not valid in the U'-n^^ rogion. It was thoroforo 
decided to utilize tost 4, primarily to investigate 
the consolidation process in this area. The use of 
•g- ton increments of loading, shov/s the pormoability 
as a linear function of the porosity up to the 6-7 
TSP range. Comparing this with Figure 17, it is noted 
that the straight lino relationship persisted only to 
the 4 TSP range, in Tests No, 2 and 3, This variance 
can only bo roasonably attributable to the influence 
of the loading increments. It is suggested that the 
instantaneous application of a major portion of the 
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total load in oach series of tests (i.e.4TSP in- 
crei.ients in Tests No< 2 and 3; 9 TSP in Test No. 4) 
causes a massive structural rearrangement v>hich takes 
place in a very short time interval. The use of 
small increments in the early portionsof both scries 
of tests produces a gradual structural rearrangement 
which takes place during a much longer time interval. 

A comparison of the relative curve shapes at 7 TSP in 
each scries of tests shows the effect of this load 
application variance. In Tests 2 and 3 the k~n curve 
in this region is definitely exponential, whereas in 
Test 4 it approaches linearity 

A pertinent question at this point is whether 
the linear relationship between k and n would continue 
through additional small load increment cycles until 
such time as a large load is applied, or whether, at 
some point, the permeability will asymtotically 
approach a condition of constancy over a largo change 
in porosity until shear failure occurs. It is 
hypothesized, from the limited test data available, 
that the latter will occur. The physical phenomenon 
involved would indicate that the gradual rearrangement 
of structure induced by small loading increments 
produces a not reduction in pore-channel effective 
area, despite the shearing of bound v/ator previously 
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discussed. In effect, the amount of shearing 
occuring, under small load increments, is not 
sufficient to retain the initial effective area and, 
thus, permeability decreases at a linear rate v/ith 
the porosity. If however wo consider the application 
of a largo load increment, the shearing of the bound 
water is of such magnitude that the effective pore- 
channel area is reduced at a much lower rate than the 
san^Dle is compressing. In effect wo have retained a 
great portion of the effective channel area despite 
a reduction in porosity. 

The extension of this phenomenon to include all 
loading, up to the point whore plastic flow will 
occur, is of course a largo jurap in the process of 
research investigation, and is subject to justifiable 
criticism, Hov/evor, it is proposed in the hope that 
future investigators in this area, may study tho 
problem in detail and by such study arrive at a 
rational analysis of the cause of tho phenomenon, 

Tho importance of detailed analysis of the k-n 
relationship is realized when one considers that tho 
accuracy of time-settlement predictions, by either 
the Terzaghi Theory or tho Schiffman Extension, 
depend upon tho porosityrpormcability relationship, 
in that tho dissipation of hydrostatic excess prossuno 
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is governed by the availability of "effective" pore 
channels for flow. Since by the hypothesis proposed 
above, this relationship is controlled not only by 
the rate of loading but also by the amount of load, 
it is apparent that further study of this problem is 
essential to a clearer analysis of the constant 
permeability ovor a finite time increment approximat- 
ion, of the Variable permeability condition, utilized 
by Sc hi ff man. 
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PART VI. 

CONCLUSIONS 

1, The assiomption of linear variation of pcrricability 

and porosity as proposed by Schmid, and utilised 
by Schiffman in developing his Extension to the 
Theory of Consolidation is valid for construction 
type loading (io. small loads applied over long 
time periods in repetitive sequence), 

2* The assumption of linear variation of permeability 
and porosity is not valid under conditions of 
large absolute load increments applied during a 
short time interval i 

3. It is hypothesized that an exponential relationship 

exists bctvifcen permeability and porosity for the 
conditions stated in 2 above. The value of the 
exponential relationship is a function not only 
of the rate of loading, but also of the magnitude 
of the load increment and it is suggested that such 
an effect is caused by the shearing of the v/atcr 
hull from the clay particles under intense 
densif ication of structure during a short time 
period. It is strongly recommended that further 
investigation bo conducted at higher load incremonts 
to tost the validity of this hypothesis'. 
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4« It is further suggested that the pori:ieahility variance, 
previeusly referred to as the "yo-yo" effect, re- 
quires that, for the finite increment approximation 
in the Variable Permeability case, developed by 
Schiffman> the increments must bo kept as small as 
possible. This limits the general usofullness of 
the approximation, since practical laboratory 
techniques preclude constant attendance. As noted 
by Schiffman (9) the degree of accuracy is a funct- 
ion of the time increments utilized. Since the 
variance of permeability is neither predictable 
nor alv/ays present, the use of the approximation 
technique for the solution of the variable 
permeability problem requires further investigation 
and analysis before a definitive conclusion can be 
drawn, 

5, No conclusion is dravifn relative to the consolidation- 
pormoability-timo relationship since fitting 
procedures for comparison of test curves and the 
theoretical curves developed by Schiffman, have 
not boon established. It is inferred from the 
permeability-porosity relationships, developed 
during this study, that the rate of loading and 
amount of load will also govern the solution of 



this problem. 
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